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The QCD corrections to the R^? moment of the invariant mass distribution in hadronic decays of the r are 
discussed. The next-to-next-to-leading order prediction is shown to be stable with respect to change of the 
renormalization scheme, provided that the contour integral expression is used. The optimized predictions are 
obtained using the principle of minimal sensitivity to select the preferred renormalization scheme. The optimized 
predictions for R}j? and R T are used in a simplified fit to the experimental data to determine the strong coupling 
constant and the parameter characterizing the nonperturbative contribution. 



Phenomenological applications of the finite or- 
der perturbative QCD predictions involve an am- 
biguity resulting from the freedom of choice of the 
renormalization scheme (RS). Formally the RS de- 
pendence is of higher order than the considered 
perturbative expression, but numerically it may 
be quite significant. In order to obtain reliable 
predictions one should therefore carefully investi- 
gate the possibilities to optimize the choice of the 
RS. Also, one should study the stability of the op- 
timized predictions by varying the RS parameters 
in some a priori acceptable range. (A condition 
on the acceptable schemes, which is applicable to 
NNLO approximants,Jaas been recently proposed 
by the present authorld.) Such an analysis is very 
important in the case of the QCD corrections to 
the hadronic r decays, for which the characteris- 
tic energy scale is m T = 1.777 GeV. In particular, 
it is interesting to what extent the RS dependent 
affects the recent results of two collaborations,l3t3'Q 
which used the hadronic r decay data to obtain a 
surprisingly precise value of a s . 

The strong interaction effects in r may now 
be studied using the QCD corrections to the R T 
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and the corrections to the R!f moments, defined 
by the relationQ 
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where T e is the electronic width of r and dT uc i/ds 
is the invariant mass distributions of the Cabbibo 



allowed hadronic decays of r. For we have: 
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where the factor Sew — 1-0194 represents cor- 
rections from electroweak interactions and Rq 1 de- 
notes the parton model predictions. For R T we 
have: 
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R T — 3SckmSew(^ + SpT + 

where S CKM = (\V ud \ 2 + \V US \ 2 ) 
contribution denotes the purely perturbative QCD 
correction, evaluated for three massless quarks. 
(We have <5£f = <5°°.) The 5 m contribution de- 
notes the correction from quark masses, which is 
practically negligible in the case of R^ 1 . For R T 
we haveH 5^ « 0.009. The S S vz contribution is 
a nonperturbative QCD, correction calculated us- 
ing the SVZ approach 
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D=4,6. 



The parameters Od in Eq.(^|) denote vacuum ex- 
pectation values of the gauge invariant operators 
of dimension D. The cd coefficients are in prin- 
ciple power series in the strong coupling constant, 
and depend on the considered moment of the in- 
variant mass distribution. 

The interest in the R!f moments comes from 
the effort to improve the accuracy of the determi- 
nation of a s from the r decay. Not unexpectedly, 
a major factor limiting the precision of the QCD 
prediction for R T is the uncertainty in the non- 
perturbative contribution. The contribution from 
the D = 4 term in the SVZ expansion for R T may 
be reliably expected to be smallB, since O4 is well 
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constrained by the sum rules phenomenology, and 
the relevant coefficient function starts at 0(a 2 s ). 
However, the D = 6 contribution to R T is not 
suppressed, and there is little information on the 
value of Og. It was therefore proposed Q to treat 
Od as free parameters which are to be extracted 
together with a s from a fit to the experimental 
data for R T and the higher moments of the in- 
variant mass distribution. Of particular interest 
here is the R} 2 moment, because similarly to R T 
the D = 4 contribution is suppressed for this mo- 
ment, and there is significant contribution from 
the D — 6 term. If in the SVZ expansion we re- 
tain only the D — 6 term, which appears to be 
a dominant source of the uncertainty in the non- 
perturbative sector, then a simplest self-consistent 
approach is to take R T and R} 2 . This is assumed 
in the following. 

A detailed discussion of the RS dependence of 
5p° l was presented elsewhere □, so in this note we 
concentrate on 5p 2 . The QCD correction S^ 2 may 
be expressed as a contour integral in the complex 
energy plane with the so called Adler function un- 
der the integral: 
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where C is a contour running clockwise from a = 
m 2 — it to a = m 2 + it away from the region 
of small \a\. In the actual calculation we assume 
that C is a circle \a\ = m 2 . The Adler function is 
defined by the relation: 

(_12^ )CT ^_ n ( / 1) (a) = 3S CK M[l + S D (-a)} (7) 
aa 

where IlL 1 ' denotes transverse part of the vector 
current correlator. The function f 12 (a/m 2 ) has 
the form: 
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The NNLO renormalization group improved per- 
turbative expansion for 5d,v may be written in 
the form: 

*D,V(-°-) = ai-a^l + nai-^ + rz^i-a)], (9) 

where a = g 2 / (Air 2 ) denotes the running coupling 
constant that satisfies the NNLO renormalization 
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Figure 1: Contour plot of 5jj as a function of the scheme 



parameters ri and C2, for 



A (3) 
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325 MeV. For technical 



reasons we use C2 — eixi on the vertical axis instead of C2- 
The boundary of the region of scheme parameters satisfying 
the Eq. fL3 is also indicated. 



group equation: 
da 



— = -- a" (1 + cia + c 2 a ). 
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In the MS scheme we have 



MS 



= 1.63982 and 
6.37101. The renormalization group co- 



efficients for ri/ =3 are b = 4.5, c\ = 16/9 and 
dp = 3863/864 w 4.471. 

By keeping the renormalization group im- 
proved expression for the Adler function under the 
integral and evaluating the contour integral nu- 
merically we obtain an essential improvement of 
the conventional perturbation expansion for 5p^, 
resumming to all orders some of the corrections 
arising from analytic xiontinuation from spacelike 
to timelike momenta. E3 

The predictions calculated in the next-to- 
next-to-leading order (NNLO) approximation de- 
pend on two RS parameters, which in principle 
may be arbitrary. Similarly as in the previous 
work 19 we parametrize the freedom of choice of 
the RS by the parameters r% and ci — the pa- 
rameter ri is determined using the RS invariant 
combination: 

92 = c 2 + r 2 - cin - r\. (11) 

For the Adler function we have pi = 5.23783. 
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To optimize the choice of the-RS we use the 
principle of minimal sensitivity Ell (PMS), which 
singles out the scheme parameters for which the 
finite order prediction is least sensitive to the 
change of RS, similarly to what we expect from 
the actual physical quantity. The dependence of 
Sp% on the scheme parameters is illustrated in fig- 
ure FD for A^L = 325 MeV. (For technical reasons 

u MS y 

we use C2 — C\T\ on the vertical axis instead of C2.) 
We choose as our optimized parameters t\ = 
and c 2 = 1.5,02 — for small values of A^L this 
point lies very close to the critical point, and even 
for large values of A^L the RS dependence in the 
vicinity of this point is very small. 

To investigate the stability of the predictions 
we use the condition proposed by the present au- 
thor, El based on the notion that natural renormal- 
ization schemes should not induce extensive can- 
cellations in the expression for the RS invariant 
P2 ■ The schemes that have the same — or smaller 
- degree of cancellation in p2 as the preferred 
scheme, satisfy in our case the inequality: 

|c2| + M+ci|n|+r?<2|p2|. (12) 

The boundary of the relevant region in the scheme 
parameter space has been denoted in the figure [l]. 
Let us note that the commonly used A/5 scheme 
lies outside of this region, but the numerical value 
of the prediction is close to the lowest value ob- 
tained within the "allowed" region. 

In figure the NNLO PMS predictions for 
Spt are shown as a function of m T /Aj^-, together 
with the minimal and maximal value obtained by 
varying the scheme parameters within the region 
determined by the condition (fl2|). We see that 
the NNLO predictions for 6%, obtained by nu- 
merically evaluating the countour integral expres- 
sion (fih, are free from potentially dangerous RS 

(3) 

instabilities even for large values of A1_L. This 

° MS 

situation is similar to that encountered for 5*°*. 
For comparison we also show the PMS predictions 
obtained in the next-to-leading order (NLO). (In 
NLO we have rf MS « -0.64.) We see that RS 
dependence of NNLO expression within the region 
defined by the condition (jl^) is smaller than the 
difference between NNLO and NLO PMS predic- 
tions. 

In order to see how the PMS optimization af- 
fects the fits to the experimental data we first test 
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Figure 2: 5i? as a function of m T /hxl obtained in the 

" MS 

PMS scheme in NNLO and NLO (upper and lower solid 
curves, respectively). The dashed lines indicate variation of 
the prediction when scheme parameters are changed within 
the region satisfying the Eq. 

the accuracy of the approximation in which one 
only retains the Oq contribution in the SVZ ex- 
pansion. To this end we make a fit of a s and Oq 
in the MS scheme, and compare the results with 
the fit performed by ALEPH0, in which the O4, 
Oq and 0% contributions have been taken into ac- 
count in the (1,0), (1,1), (1,21 and (1,3) moments. 
If we take, following ALEPHB, R t = 3.645 ±0.024 
and D\ 2 = R^/Rf = 0.0570 ± 0.0013, we ob- 
tain from the fit in the MS scheme a s (M§) = 
0.1209 ± 0.0013 and O e = -0.0010 ± 0.0012. This 
appears to be remarkably close to the values 0.121 
and -0.0016 obtained in the full fit by ALEPH. 
This gives us confidence that the 11 Oq approxima- 
tion" captures the essential features of QCD cor- 
rections in r decays. 

Performing the same fit, but using now the 
NNLO PMS predictions, we obtain a s (M§) = 
0.1198 and 6 = -0.0011, i. e. the value of the 
condensate practically does not change, but the 
value of the strong coupling constant is reduced by 
about one standard deviation of the experimental 
error. It is also of some interest to compare these 
results with the NLO PMS fit — we then obtain 
a s (M§) = 0.1221. Taking the difference of the 
NNLO and NLO PMS fits is perhaps the best way 
of estimating the accuracy of the perturbative pre- 
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diction. We see that thus obtained uncertainty of 
a s is of the order 0.0022, i. e. it is quite large, 
compared for example to the experimental uncer- 
tainty. 

We may now use the optimized predictions to 
obtain a s from the more up to date experimental 
data on r decays. We take R T — 3.613 ± 0.032, 
which is a weighted average of three possible de- 
terminations involving B e = 0.1790 ± 0.0017, 

= 0.1744 ± 0.0023 and t t = (292.0 ±2.1) x 
10~ 15 sec (values according tp-.the 1995 update 
of the Particle Data Group E3). We also take 
Dl 2 = 0.0561 ± 0.000J6, which is a weighted av- 
erage of the ALEPH and CLEO □ determina- 
tions. With these numbers we obtain a s (M|) = 
0.1177 ± 0.0017 K(m2) = 0.321 ± 0.014) and 
O e = -0.0020 ±0.0005. 

Concluding, the perturbative QCD correction 
to the R 12 moment of the invariant mass distri- 
bution in hadronic tau decays was found to be 
stable with respect to change of the renormaliza- 
tion scheme, despite the low energy scale, pro- 
vided that the contour integral expression is used, 
which resums some of the corrections to all orders. 
However, the difference between predictions in the 
conventionally used MS scheme and the PMS 
scheme was found to be phenomenologically signif- 
icant. Also, the difference between the NNLO and 
NLO predictions in the PMS scheme was found 
to be significant, indicating perhaps that the un- 
certainty in the perturbative prediction is larger 
than previously expected. Finally, the optimized 
predictions have been used to obtain a realistic fit 
for a s from the experimental data. 
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